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ABSTRACT 

There is currently debate over whether the dust content of planetary systems is stochastically 
regenerated or originates in planetesimal belts evolving in quasi-steady state. In this paper a simple 
model for the steady state evolution of debris disks due to collisions is developed and confronted with 
the properties of the emerging population of 7 sun-like stars that have hot dust at < 10 AU. The model 
shows that there is a maximum possible disk mass at a given age, since more massive primordial disks 
process their mass faster. The corresponding maximum dust luminosity is / max = 0.16 x 10~ 3 r 7 / 3 £~ g 1 e , 
where r is disk radius in AU and i a ge is system age in Myr. The majority (4/7) of the hot disks 
exceed this limit by a factor 3> 1000 and so cannot be the products of massive asteroid belts, rather 
the following systems must be undergoing transient events characterized by an unusually high dust 
content near the star: r/ Corvi, HD69830, HD72905 and BD+20307. It is also shown that the hot 
dust cannot originate in a recent collision in an asteroid belt, since there is also a maximum rate at 
which collisions of sufficient magnitude to reproduce a given dust luminosity can occur in a disk of 
a given age. For the 4 transient disks, there is at best a 1:10 chance of witnessing such an event 
compared with 2% of stars showing this phenomenon. Further it is shown that the planetesimal belt 
feeding the dust in these systems must be located further from the star than the dust, typically at 
3> 2 AU. Other notable properties of the 4 hot dust systems are: two also have a planetesimal belt 
at > 10 AU (17 Corvi and HD72905); one has 3 Neptune mass planets at < 1 AU (HD69830); all 
exhibit strong silicate features in the mid-IR. We consider the most likely origin for the dust in these 
systems to be a dynamical instability which scattered planetesimals inwards from a more distant 
planetesimal belt in an event akin to the Late Heavy Bombardment in our own system, the dust being 
released from such planetesimals in collisions and possibly also sublimation. Further detailed study of 
the planet, planetesimal and dust populations in these rare objects has the potential to uncover the 
chaotic evolutionary history of these systems and to shed light on the history of the solar system. 
Subject headings: circumstellar matter — planetary systems: formation 



1. INTRODUCTION 

Planetesimal belts appear to be a common feature of 
planetary systems. There are two main belts in the solar 
system: the asteroid belt and the Kuiper belt. These 
belts inhabit the regions of the solar system where plan- 
etesimal orbits can remain stable over the 4.5 Gyr age 
of our system (Lecar et al. 2001). The larger planetes- 
imals in the belts are continually grinding down feeding 
the smaller bodies in a process known as a collisional 
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cascade which is slowly eroding the belts (Bottke et al. 
2005). The smallest dust in the asteroid belt is acted on 
by radiation forces; P-R drag makes the dust spiral in to- 
ward the Sun making a disk known as the zodiacal cloud 
that the Earth sits in the middle of (Leinert & Grim 
1990). A dust cloud is also predicted to arise from colli- 
sions amongst Kuiper belt objects (Liou & Zook 1999), 
although our information on this population is sparse 
(Landgraf et al. 2002) because its emission is masked 
by the zodiacal emission (Backman, Dasgupta & Stencel 
1995) and few dust grains make it into the inner solar 
system (Moro-Martm & Malhotra 2003). 
Many extrasolar systems also have such planetesimal 
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belts, known as debris disks. These have been detected 
from their dust content (Aumann et al. 1984) from which 
it has been inferred that larger planetesimals must exist 
to replenish the dust disks because of the short lifetime 
of this dust (Backman & Paresce 1993). The collisional 
cascade scenario is supported by modeling of the emis- 
sion spectrum of the dust which shows a size distribution 
similar to that expected for dust coming from a colli- 
sional cascade (Wyatt & Dent 2002, hereafter WD02). 
However, the issue of how these disks evolve has recently 
come under close scrutiny. 

From a theoretical point view, Dominik & Decin (2003; 
hereafter DD03) showed that if P-R drag is not important 
then a planctcsimal belt evolving in quasi-steady state 
would lose mass due to collisional grinding down giving a 
disk mass (and dust luminosity) that falls off oc This 
is in broad agreement with the observed properties of 
debris disks: the mean dust luminosity at a given age falls 
off oc t -18 (Spangler et al. 2001); the mass inferred from 
detection statistics falls off oc t~ - 5 (Greaves & Wyatt 
2003), while the mass of the detected disks falls off oc t~ x 
(Najita & Williams 2005): the upper limit in luminosity 
of the detected disks also falls off oc i -1 (Rieke et al. 
2005). While these trends can be viewed as a success of 
the steady-state model, it has yet to be proved that a 
steady state evolution model fits the data in more than 
just general terms (Meyer et al. 2006). Several puzzling 
observations also remain to be explained. 

Decin et al. (2003) noted that the maximum frac- 
tional luminosity of debris disks remains constant at 
/ = Li-c/Li, rts 10~ 3 up to the oldest stars, where L a 
and L* are the disk and stellar luminosities respectively 
(see also Table Al for definitions of the parameters used 
in the text), and this was explained by DD03 as a con- 
sequence of delayed stirring. A delay in the ignition of 
a collisional cascade is expected if it is the formation of 
Pluto-sized objects which trigger the cascade, since such 
massive bodies take longer, up to several Gyr, to form 
further from the star (Kenyon & Bromley 2002). How- 
ever, that interpretation predicts that the radius of the 
belts should increase with stellar age, and this is not 
observed (Najita & Williams 2005). There is also re- 
cent evidence that the dust content of some systems is 
transient. The discovery of a population of dust grains 
around Vega in the process of removal by radiation pres- 
sure indicates that this system cannot have remained in 
steady state for the full 350 Myr age of the star (Su et 
al. 2005). Rieke et al. (2005) used their statistics on A 
stars, which showed a wide variety of properties among 
the debris disks, to suggest that much of the dust we 
see is produced episodically in collisions between large 
planetesimals. There is also an emerging population of 
debris disks detected around sun-like stars with dust at 
a few AU (Gaidos 1999; Beichman et al. 2005; Song et 
al. 2005; Smith, Wyatt & Dent in prep.). There is de- 
bate over whether these are atypically massive asteroid 
belts or the consequence of a rare transient event (e.g., 
Beichman et al. 2005). 

A stochastic element to the evolution of debris disks 
would fit with our understanding of the evolution of the 
dust content of the inner solar system. This is believed 
to have been significantly enhanced for timescales of a 
few Myr following collisions between objects ^100 km 
in size in the asteroid belt (Nesvorny et al. 2003; Far- 



ley et al. 2006). However, it is not known whether the 
aftermath of individual collisions would be detectable in 
a debris disk, or indeed whether such events would hap- 
pen frequently enough to explain the statistics (WD02; 
Telesco et al. 2005). Such events have a dramatic effect 
on the amount of dust in the solar system because there 
is relatively little around during the quiescent periods. 
Planctcsimal belts of equivalent mass to those in the so- 
lar system would not have been detected in the current 
debris disk surveys. However, there is evidence to sug- 
gest that both belts were ~ 200 times more massive in 
the past (e.g., Stern 1996; Bottke et al. 2005). Periods 
analogous to the heavy bombardment experienced in the 
solar system up to ~ 700 Myr after its formation have 
also been invoked to explain the fact that debris disks are 
most often detected around stars < 400 Myr old (Habing 
et al. 1999). 

In the light of this controversy we revisit a simple an- 
alytical model for the steady state collisional evolution 
of planetesimal belts which was originally explored in 
DD03. The model we derive for that evolution is given 
in §2, and differs in a subtle but important way from that 
of DD03, since it affects the dust production as a func- 
tion of collision velocity. This model shows that there 
is a maximum possible disk mass (and dust luminosity) 
at any given age. In §3 confrontation with the few hot 
planctcsimal belts discovered recently shows that the ma- 
jority of these cannot be explained as massive asteroid 
belts, rather these must be systems undergoing a tran- 
sient event. The possibility that these are caused by 
a recent collision within a planctcsimal belt is also dis- 
cussed, as is the possibility that the dust originates in 
a planetesimal belt in the terrestrial planet region. The 
implications of these results are discussed in §4. Appli- 
cation of the model to the statistics of detected debris 
disks will be considered in a later paper (Wyatt et al. , in 
prep.). 

2. ANALYTICAL COLLISIONAL EVOLUTION MODEL 

In this section a simple analytical model is developed 
for the evolution of a planetesimal belt due to collisions 
amongst its members. The parameters used in this model 
are summarized in the table Al which also gives the units 
assumed for these parameters throughout the paper. 

2.1. The planetesimal belt size distribution 

The planetesimal belt is assumed to be in collisional 
equilibrium with a size distribution defined by: 

n(D) = KD 2 - 3q , (1) 

where q = 11/6 in an infinite collisional cascade 
(Dohnanyi 1969) and the scaling parameter K is called f a 
by DD03. That distribution is assumed to hold from the 
largest planetesimal in the disk, of diameter D c , down to 
the size below which particles are blown out by radiation 
pressure as soon as they arc created, £>bi- If we assume 
that q is in the range 5/3 to 2 then most of the mass is in 
the largest planetesimals while the cross-sectional area is 
in the smallest particles such that: 

a tot =3.5x W-^K^q-^y^lO^Du) 5 -^ (2) 

M tot = 8.8 x 10- 17 Kp(6 - Zq)- 1 !^- 3 *, (3) 
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where spherical particles of density p have been assumed 
and Mtot is in if the units of table Al are used for 
the other parameters. 

The planetesimal belt is assumed to be at a radius r, 
and to have a width dr (in AU). One of the observable 
properties of a planetesimal belt is its fractional lumi- 
nosity, / = Ljr/L*, i.e., the infrared luminosity from the 
disk divided by the stellar luminosity. Assuming that 
the grains act like black bodies and so absorb all the 
radiation they intercept we can write: 

/ = <7tot/(47IT 2 ). (5) 

In other words, in this model <7tot, A^tot and / are all 
proportional to each other and just one is needed to de- 
fine the scaling factor K in equation (1). Assuming the 
particles act like black bodies also allows us to derive the 
following relation: 



A,i = 0.8(L*/M*)(2700/p), 



(6) 



where D D i is in fim, L± and M* are in solar units, and p 
is in kg m~ 3 . 

Relaxing the black body assumption is easily achieved 
(e.g., WD02). However, this would result in relatively 
small changes in the way / scales with M to t, and so for 
its heuristic simplicity we keep this assumption through- 
out this paper. Probably the most important simplifi- 
cation within this model is that of the continuous size 
distribution. For example, we know that the cut-off in 
the size distribution at -Dbi would cause a wave in the 
size distribution at sizes just larger than this (Thebault, 
Augereau & Beust 2003), that large quantities of blow- 
out grains can also affect the distribution of small size 
particles (Krivov, Mann, & Krivova 2000), and that the 
dependence of planetesimal strength on size can result in 
q =/= 11/6 as well as a wave in the distribution at large 
sizes (Durda et al. 1998; O'Brien & Greenberg 2003). 
Also, since the largest planetesimals would not be in col- 
lisional equilibrium at the start of the evolution, their 
initial distribution may not be the same as that of a col- 
lisional cascade, although distributions with q w 11/6 
have been reported from planet formation models (e.g., 
Stern & Colwell 1997; Davis & Farinella 1997; Kenyon 
& Luu 1999) meaning this is a reasonable starting as- 
sumption. Despite these simplifications, we believe this 
model is adequate to explore to first order the evolution 
of planetesimal belts which can later be studied in more 
depth. 

2.2. Collisional evolution 

In a collisional cascade material in a bin with a given 
size range D to D + dD is replaced by fragments from 
the destruction of larger objects at the same rate that 
it is destroyed in collisions with other members of the 
cascade. The long-timescale evolution is thus determined 
by the removal of mass from the top end of the cascade. 
In this model the scaling factor K (and so the total mass 
and fractional luminosity etc) decreases as the number of 
planetesimals of size D c decreases. The loss rate of such 
planetesimals is determined by their collisional lifetime, 
which in the terminology of WD02 is given by: 

t c = Vr 3 /M*(rdr/ CTt „ t )[2///(e, /)]// cc , (7) 

where maintaining the units used previously gives t c in 
years, I is the mean inclination of the particles' orbits 



(which determines the torus height), f(e,I) is the ratio 
of the relative velocity of collisions to the Kcplcrian ve- 
locity (= Wroi/^k, also called v by DD03), and / cc is the 
fraction of the total cross-sectional area in the belt which 
is seen by planetesimals of size D c as potentially causing 
a catastrophic collision. 

From hereon we will use the assumption that f(e, I) — 
\/l-25e 2 + I 2 , where e is the mean eccentricity of the par- 
ticles, which is valid for Rayleigh distributions of e and 
I (Lissauer & Stewart 1993; Wetherill & Stewart 1993). 
An expression for f cc was given in WD02, however, here 
we will ignore the gravitational focussing effect, which is 
important in the accumulation phase but not during the 
destruction phase of a planetesimal belt (see §3.2), and 
so derive an expression that is the same as that given in 
Wyatt et al. (1999): 

/cc = (10- 9 D hl /D c )^- 5 G(q,X c ), (8) 

where X c = D cc /D c , D cc is the smallest planetesimal 
that has enough energy to catastrophically destroy a 
planetesimal of size D c (which is called e in DD03), and: 

G(q,X c ) = {(Xt 3q - 1) + (6g - 10)(3g - A)~ 1 (Xt 3q - 1) 

+ (3q-5)(3q-3)- 1 (Xt 3q ^l)}- (9) 

The factor X c can be worked out from the dispersal 
threshold, defined as the specific incident energy re- 
quired to catastrophically destroy a particle such that 
(WD02): 

X c = (2Qyvl l ) 1 / 3 , (10) 

= 1.3 x lO-^rM-Vte, iy 2 } 1 ' 3 , (11) 

where Q* p is in J kg" 1 (called S in DD03 2 ) 

Combining the above equations gives for the collisional 
lifetime of the planetesimals of size D c : 

5-3q 

= ( ___ .. ] ( - - ! ( - ) (12) 



.(13) 



3.8pr 25 drD, 



tot 



G(q,X c ) J V D, 

(12 g -20)[l + 1.25(e/7) 2 ]-°- 5 
(18-9q)G(q,X c ) 



Assuming that collisions are the only cause of mass loss 
in the belt, the evolution of the disk mass M tot (t) (or 
cquivalently of K, otot, or /) can be worked out by solv- 
ing dM tot /dt = -M tot /t c to give: 

M tot (t) = M tot (0)/[1 + i/*c(0)], (14) 

where M to t(0) is the initial disk mass and i c (0) is the 
collisional lifetime at that initial epoch; this solution is 
valid as long as mass is the only parameter of the plan- 
etesimal belt that changes with time. This results in a 
disk mass which is constant at M to t(0) for t <C t c (0), 
but which falls off oc l/t for t S> t c (0) (as noted, e.g., in 
DD03). 

However, another interesting property of this evolution 
is that, since the expression for t c (0) includes a depen- 
dence on Mtot(0), the disk mass at late times is indepen- 
dent of initial disk mass. This is because more massive 

2 Equation 25 in DD03 differs from our equation (10) because 
we define Q£> to be the specific incident kinetic energy so that 
0.5M2f 2 ol = MiQq whereas DD03 define S to be the specific bind- 
ing energy of the two objects (giving their equation 24). In the limit 
of S <?C « 2 el /8 the two equations are the same, since l c < 1. 
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Fig. 1. — The dependence of (left) G(11/6,X C ) and (right) X c on planetesimal eccentricity (e) for planctcsimals of different strengths 
(Qd) an d at different distances from the star (r). 



disks process their mass faster. This means that for any 
given age, t agc , there is a maximum disk mass M max 
(and also infrared luminosity, / max ) that can remain due 
to collisional processing: 



M max = 



3.8 x 1(T 6 pr 3b {dr/r)D, 



age 



21-0.5 



(12g- 20) [1 + 1.25(e/J) 
(18- 9q)G(q,X, 

10- 6 r 1 - 5 (d r/r) \ ^10" 9 £> bl 



ago 



D c 



2[l + 1.25(e/J) 2 
G(q,X c ) 



(15) 



(16) 



In this model, the present day disk mass (or luminosity) 
is expected to be equal to this " maximum" disk mass (or 
luminosity) for disks in which the largest planetesimals 
arc in collisional equilibrium. This corresponds to disks 
around stars that are older than the collisional lifetime 
of those planetesimals given in equation (13). 

For example, with the further assumptions that q = 
11/6, e w I, and p = 2700 kg m~ 3 , wc find: 

M max = 0.009r 3 - 5 (dr/r)^ c M-°- 5 C o /G(ll/6,X c )(17) 



age/ 

/ max = 0.004r 1 - 5 (rfr/r)i? c - 5 L; - 5 t ag 1 c /G'(ll/6,X c I18) 



where M max is in Mq, r in AU, D c in km, i age in Myr, 
and G(U/6,X C ) = X~ 5 + 0.67X- 15 + 0.2X~ 2 5 - 1.87, 

with X c = 10- 3 (rQ^/e 2 ) 1/3 (<9d is in J kg" 1 ). 

Plots of G(11/6,X C ) and X c for typical planetes- 
imal belts are shown in Fig. 1. However, for 
many disks the approximation that X c <C 1 is 
valid, and so G(11/6,X C ) w 0.2X~ 2 - 5 = 6.3 x 

10 6 r- 5 / 6 Q^- 5/6 e 5 / 3 M* 5/6 , giving: 



Mr, 



: = 1.4 x 10- 9 r 13 / 3 (dr/r)D c Q* D 



5/6 



e- 5/3 M: 4/3 C. (19) 
/max -0.58 x 10-V/ 3 (dr/r)^ c °- 5 Q^ 5/6 x 

e-^M-^L-O-^. (20) 



2.3. Comparison with DD03 



Since DD03 produced a very similar analytical model, 
our results were compared with those of DD03. The 
results of disk evolution for a planetesimal belt close 
to their nominal model were computed using the pa- 
rameters: r = 43 AU, dr = 15 AU, D c = 2 km, 
p = 2700 kg m- 3 , f(e,I) = 0.1, e/I = 1, Q* = 200 J 
kg" 1 , M tot (0) = 10M®, AO star (for which U = 54L Q , 
M+ = 2.9M Q , Z?bi = 15 /an). Each of the parameters 
A^tot(O), r, f(e, I), D c and spectral type were also varied 
to make the plots shown in Fig. 2 which are equivalent 
to Figs lb-lf of DD03. 

The results are very similar in most regards: more 
massive disks start out with higher /, but the turnover 
from constant to 1/t evolution is later for lower mass 
disks meaning that at late times all disks converge to the 
same maximum value (Fig. 2a); putting the same mass 
at larger distances reduces the initial dust luminosity /, 
but the resulting lower surface density and longer orbital 
timcscales there combine to make the turnover happen 
later which means that at late times more distant belts 
arc more massive (Fig. 2b); putting the same mass into 
larger planctcsimals reduces the cross-sectional area of 
dust (equation 4) and so the initial dust luminosity /, 
but increases the collisional lifetime of those planetesi- 
mals (equation 13) which means that at late times belts 
with larger planctcsimals retain their mass for longer 
(Fig. 2d); later spectral types have higher starting dust 
luminosities because the cascade extends down to smaller 
sizes (equation 6), and the longer orbital times mean that 
they keep their mass for longer (Fig. 2e). 

Where the models differ is in the exact way M tot is 
used to get / and t c , and in the way the evolution is af- 
fected by changing w ro i/wk (Fig. 2c). This is because the 
models make different assumptions. Here we assume that 
the size distribution is continuous between D c and Dbi, 
whereas in DD03 the large planetesimals feeding the cas- 
cade are seen as separate from the cascade. This means 
that for us M tot gives a direct estimate of K (equation 
3) and so the amount of dust /, while for DD03 they 
equate the mass flow through the cascade with the mass 
input from the break-up of planetesimals meaning that 
while their scaling parameter is proportional to M tot (as 
is ours), it also includes a dependence on the parameter 
we call X c which affects the mass flow rate in the cas- 
cade. This explains all of the differences: the details of 
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Fig. 2. — The collisional evolution of a planetesimal belt with parameters similar to the nominal model of DD03 [r = 43 AU, dr = 15 
AU, D c = 2 km, p = 2700 kg m~ 3 , f(e,I) = 0.1, e/I = 1, Q* = 200 J kg" 1 , M tot (0) = 10M e , AO star] showing the effect of changing: 
(top left) starting disk mass M tot (0), (top right) disk radius r, (middle left) collision velocity v Ie i/v\^, (middle right) maximum 
planetesimal size D c , and (bottom left) stellar spectral type. These plots can be directly compared to figs, lb-f of DD03. 



the scaling explain the slightly different initial / values 
in all the figures, and the fact that for us planetesimals of 
size D c are destroyed by planetesimals down to size X C D C 
means that our collisional lifetimes are always shorter 
than those in DD03, since they assume that planetes- 
imals only collide with same size planetesimals. For us 
changing u r ci/wk does not affect the initial / parameter as 
described above, but it does affect the collisional lifetime 
of the largest planetesimals which can survive longer if 
Wrci/^k is reduced (since this means that fewer planetes- 
imals in the cascade cause destruction on impact). The 
opposite is the case for the DD03 model: changing v Te \/vk 
does not affect the collisional lifetime of the largest plan- 



etesimals, since they only collide with each other, but a 
lower collision velocity does increase the initial dust lu- 
minosity because the cascade must have more mass in it 
to result in a mass flow rate sufficient to remove mass in- 
troduced by the large planetesimals. While the difference 
is subtle, it is important, since f re i/uk may be important 
in determining the presence of dust at late times (DD03; 
section 3). 

On the face of it, it seems that our model provides a 
more accurate description of the disk. The reason is that 
in a collisional cascade the mass flow does not need to be 
taken into account, since it results in the q = 11/6 size 
distribution (Tanaka et al. 1996). In other words the 
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dependence of the scaling of the cascade with X c found 
by DD03 should have been removed if the largest plan- 
etesimals had been allowed to collide with smaller plan- 
etesimals (since increasing X c would have both restricted 
mass flow within the cascade and slowed down the mass 
input from the destruction of large planetesimals) . How- 
ever, it is also true that the q = 11/6 distribution only 
applies in an infinite cascade, and since both models have 
truncated the size distribution at D Cl this would affect 
the evolution. Also, the effect of the variation of Q£, with 
D on the size distribution and its evolution are not yet 
clear, and neither is the evolution of the size distribution 
while the collisional cascade is being set up. These is- 
sues will be discussed only briefly in this paper, in which 
the simple evolution model described above is applied to 
some of the latest observational results on debris disks. 

3. APPLICATION TO RARE SYSTEMS WITH HOT DUST 

Very few main sequence stars exhibit hot dust within 
~ 10 AU, i.e., in the region where we expect planets may 
have formed. Four surveys have searched for hot dust 
around sun- like stars (main sequence F, G or K stars) by 
looking for a 25 /im flux in excess of photospheric levels 
using IRAS (Gaidos 1999), ISO (Laureijs et al. 2002) and 
Spitzer (Hines et al. 2006; Bryden et al. 2006). All con- 
cluded that only 2 ± 2% of these stars have hot dust with 
infrared luminosities / = L^/L^ > 10~ 4 , finding a total 
of 3 candidates. Other hot dust candidates exist in the 
literature, however some IRAS excess fluxes have turned 
out to arise from chance alignments with background ob- 
jects (e.g., Lisse et al. 2002), including the candidate 
HD128400 from the hot dust survey of Gaidos (1999) 
(Zuckcrman, priv. comm.). Thus confirmation of the 
presence of dust centred on the star using ground- and 
space-based mid-IR imaging is vitally important (Smith, 
Wyatt & Dent, in prep.). The tally of confirmed hot dust 
sources now stands at seven, and these are summarized 
in Table 1 which also gives the estimated radial location 
of the dust based on fitting of the spectral energy dis- 
tribution of the excess emission; for all stars the dust is 
predicted to lie at < 10 AU. 

While the frequency of the presence of such emission 
is low, there is as yet no adequate explanation for its ori- 
gin and why it occurs in so few systems. Analogy with 
the solar system suggests that these are systems in which 
we are witnessing the collisional grinding down of atyp- 
ically massive asteroid belts. However, other scenarios 
have also been proposed in which the dust is transient, 
having been produced in some stochastic process. Such 
a process could be a recent collision between two massive 
protoplanets in an asteroid belt (Song et al. 2005), the 
sublimation of one supercomet (Bcichman et al. 2005), 
or the sublimation of a swarm of comets, possibly scat- 
tered in from several tens of AU in an episode analogous 
to the period of Late Heavy Bombardment in the solar 
system (Gomes et al. 2005). 

3.1. Are these massive asteroid belts? 

Here we consider the possibility that these are atypi- 
cally massive asteroid belts, and show that for the major- 
ity of the known systems this is unlikely to be the case. 
The reason is that given in §2.2, which is that more mas- 
sive asteroid belts are not necessarily more dusty at late 
times, and there is a maximum dust luminosity we can 



expect for a belt of a given age, given its radial loca- 
tion (equations 15-20). To arrive at a rough estimate of 
the maximum possible / max we assume the following pa- 
rameters: the largest possible planetesimal is D c = 2000 
km, since this is above the largest members of the aster- 
oid and Kuiper belts, and fits with the expectation that 
planetesimal growth is halted once the largest planetes- 
imals reach this size due to the resulting gravitational 
perturbations (Kenyon & Bromley 2002); belt width is 
air = 0.5r; planetesimal strength is = 200 J kg -1 , the 
canonical value used in DD03, although gravity strength- 
ening can give rise to higher values for planetesimals 
larger than ~ 1 km (see §3.2); planetesimal eccentricity 
is e = 0.05, typical for planetesimal belts like the asteroid 
belt that are undergoing a collisional cascade, and close 
to that expected from stirring by 2000 km planetesimals 
within such belt. 3 Substituting in these nominal val- 
ues into equation (20) and approximating M+ = L+ = 1 
gives: 

/ max = 0.16 x lO-V/^" 1 ,. (21) 

Plots analogous to those in Fig. 2 are presented in Fig. 3 
which shows the evolution for a planetesimal belt with 
the nominal parameters described above (and with a 
nominal starting mass of M to t(0) = l-M©) along with the 
consequence for the evolution of changing any of those 
parameters. Note that it is most appropriate to refer to 
Fig. 3, rather than Fig. 2, when considering the evolution 
of planetesimal belts close to sun-like stars. 

The value of / max is quoted in Table 1 under the as- 
sumption that the planetesimal belt has the same age 
as the star. The quoted value for each star is that from 
equation (18) for its spectral type, but is within a fac- 
tor of three of that given in equation (21), indicating 
that this equation may be readily applied to observed 
belts in the future. The four oldest systems (BD+20307, 
HD72905, 7] Corvi and HD69830) have / obs > 10 3 / max . 
We show in §3.2 that even with a change in parameters 
it is not possible to devise asteroid belts in these systems 
that could survive to the age of the stars giving rise to 
the observed dust luminosities. Thus we conclude that 
this period of high dust luminosity started relatively re- 
cently The timescale over which a belt can last above 
a given luminosity, / obB , is i a g e /max//obs, since collisions 
would grind a belt down to this level on such a timescale. 
This implies that belts this luminous only last between a 
few thousand years (BD+20307 and HD72905) and a few 
Myr (77 Corvi and HD69830). However, the true duration 
of this level of dust luminosity depends on the details of 
the process causing it, and moreover there is still up to 
two orders of magnitude uncertainty in / max (see §3.2). 
Thus this calculation should not yet be used to infer from 
the <~ 2% of systems with hot dust that, e.g., every sun- 
like star must undergo 10-1000 such events in its lifetime 
(or fewer systems must undergo even more events). For 
now the conclusion is that these systems cannot be plan- 
etesimal belts that have been evolving in a collisional 
cascade for the full age of the star. 

This leaves open the possibility that the collisional cas- 
cade in these systems was initiated much more recently, 

3 Equating the velocity dispersion in the belt with the es- 
cape velocity of a planetesimal of size D c gives e 2.6 X 
10 -7 p 0.5 r 0.5 M -0.5 D ^ 
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TABLE 1 

Main sequence sun-like (F, G and K) stars in the literature with evidence for hot dust at < 10 AU. 



Star name 


Sp. Type 


Age, Myr 


Radius, AU 


/obs — L\ T /L+ 


/max 


Transient? 


Reference 


HD98800 C 


K4/5V 


~ 10 


2.2 


220 x 10~ 3 


270 x 10~ 6 


Not req 


Low ct al. (2005) 


HD113766 ac 


F3V 


16 


3 


2.1 x 10~ 3 


45 x 10~ 6 


Not req 


Chen et al. (2005) 


HD12039 


G3/5V 


30 


4-6 


0.1 x 10~ 3 


200 x 10~ 6 


Not req 


Hines ct al. (2006) 


BD+20307 a 


GOV 


300 


1 


40 x 10~ 3 


0.36 x 10~ 6 


Yes 


Song ct al. (2005) 


HD72905 a 


G1.5V 


400 


0.23 b 


0.1 x 10~ 3 


0.011 x 10~ 6 


Yes 


Beichman ct al. (2006a) 


r] Corvi a 


F2V 


1000 


l-2 b 


0.5 x 10~ 3 


0.15 x 10~ 6 


Yes 


Wyatt et al. (2005) 


HD69830 a 


KOV 


2000 


1 


0.2 x 10~ 3 


0.13 x 10" 6 


Yes 


Beichman ct al. (2005) 



a infrarcd silicate feature 

b also has cool dust component at > 10 AU 
c binary star 



perhaps because a long timescale was required to form 
the 2000-3000 km sized planetesimals necessary to stir 
the planetesimal belt and cause the switch from accre- 
tion to collisional cascade (Kenyon & Bromley 2004). 
However, we consider this to be unlikely, because the 
timescale for the formation of objects of this size at 1 
AU from a solar mass star was given in Kenyon & Brom- 
ley (2004) to be ~ 0.6dr/M tot Myr, where M tot is the 
mass of material in an annulus of width dr, just as in the 
rest of the paper. This means that the cascade can only 
be delayed for 100-1000 Myr at 1 AU for planetesimal 
belts of very low mass, which would also be expected to 
have low dust luminosities when the cascade was even- 
tually ignited. For example, a delay of > 500 Myr would 
require < 0.6 x 10~ 3 M® in the annulus at 1 AU of 0.5 
AU width, a mass which corresponds to a fractional lu- 
minosity of < 5 x 10~ 6 (equations 4 and 5 with p = 2700 
kg m~ 3 and q = 11/6), much lower than that observed in 
all systems. One can also consider the same argument in 
the following way: the observed luminosity / OD s implies 
a planetesimal belt mass which current planet formation 
theories indicate would result in the growth of 2000 km 
planetesimals which would ignite a collisional cascade on 
a timescale of 3 x 10~ 3 (dr / V) / / bs Myr if this was placed 
at 1 AU from a solar mass star. The conclusion at the 
end of the last paragraph also considers the collisional 
cascade to evolve in quasi-steady state, and it is possi- 
ble that collisions between large members of the cascade 
may have recently introduced large quantitites of small 
dust; that possibility is discussed in §3.3. 

For the three youngest systems the conclusions are less 
clear. The dust luminosities of HD12039 and HD113766 
are, respectively, close to and fifty times higher than the 
maximum allowed value for collisionally evolved plan- 
etesimal belts. However, given the uncertainties in the 
parameters in the model (described in §3.2), we conclude 
that it is not possible to say that these could not be 
massive asteroid belts. The main reason that firm con- 
clusions cannot be drawn is the large radial location of 
the dust at > 2 AU. The strong dependence of / max 
on r means that it is easiest to constrain the nature of 
belts within a few AU which evolve very rapidly. For 
the youngest system (HD98800), while its dust luminos- 
ity lies a factor of 800 above the maximum for the age 
of the star, we do not infer that this must be transient, 
since the high dust luminosity and low age imply that 
this system is in a transitional phase and the collisional 
cascade in this debris disk is likely to have only recently 



been ignited. Rather we note that this model implies 
that due to collisional processing this debris disk can- 
not maintain this level of dust emission beyond the next 
<~ 10,000 years (albeit with an additional two orders of 
magnitude uncertainty, §3.2). 

3.2. Possible caveats 

Given the large number of assumptions that went into 
the estimate for / max , it is worth pointing out that this 
model is in excellent agreement with the properties of the 
asteroid belt in the solar system, since for a 4500 Myr 
belt at 3 AU the model predicts M max = 0.4 x 10~ 3 M e , 
which is close to the inferred mass of the asteroid belt of 
0.6 x 10" 3 M e (Krasinsky et al. 2002). The model also 
predicts / max = 5 x 10~ 7 , which is consistent with the 
estimate for the zodiacal cloud of L^/L* = 0.8 x 10~ 7 
(Backman & Paresce 1993). 4 It is also necessary to ex- 
plore if there is any way in which the parameters of the 
model could be relaxed to increase / ma x and so change 
the conclusions about the transience of the hot dust sys- 
tems. Equation (20) indicates one way in which / max 
could be increased, which is by either reducing the ec- 
centricities of the planetesimals, e, or increasing their 
strength, Qp, both of which could increase X c and so 
decrease the rate at which mass is lost from the cascade 
(e.g., fig. 3). The other way is to change the size distri- 
bution so that a given disk mass results in a significantly 
larger dust luminosity, e.g., by increasing q. 

In fact Benz & Asphaug (1999) found a value of 
that is higher than 2 x 10 5 J kg -1 for planetesimals 
as large as 2000 km for both ice and basalt composi- 
tions. This would result in an increase in / max by a 
factor of ~ 170 (e.g., Fig. 3). However, such a high value 
of is possible only due to gravity strengthening of 
large planetesimals, and the dependence in this regime 
of oc D 13 (Benz & Asphaug 1999) would result in an 
equilibrium size distribution with q g w 1.68, since when 
Qjf, oc D s then q = (11 + s)/(6 + s) (O'Brien & Green- 
berg 2003). If such a distribution was to hold down to 
the smallest dust grains the net result would be a de- 
crease in / max by ~ 200. This is not the case, how- 

4 In planetesimal belts as tenuous as the asteroid belt, the effect 
of P-R drag is important (Wyatt 2005) meaning that the cross- 
sectional area of dust in the zodiacal cloud is dominated by ~ 
100 fim sized grains rather than grains of size D\>\ as assumed in 
the simple model of §2.1. Taking this into account would reduce 
the fractional luminosity predicted by the model by an order of 
magnitude. 
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Fig. 3. — The collisional evolution of a planetesimal belt at r = 1 AU around a sun-like star (L* = M* = 
assuming that belt can be described by the parameters used in §3.1 (i.e., dr/r = 0.5, D c = 2000 km, p ■■ 

= 200 J kg - 1 ). All panels show dust luminosity / = L; r /L* as a function of time, and the evolution with the above nominal parameters 
is shown with a solid line. The different panels show the effect of changing the following parameters: (top left) starting disk mass Mt o t(0), 
(top right) disk radius r, (middle left) planetesimal eccentricity e, (middle right) maximum planetesimal size D c , (bottom left) 
planetesimal strength Q^, and (bottom right) stellar spectral type. 

ever, since objects in the size range D < D t — 0.15 km 
arc in the strength scaled regime where oc D 
leading to a size distribution with q s = 1.89 in this 
range. According to O'Brien & Greenberg (2003) the 
size distribution of a collisional cascade with a realistic 
prescription should have two components (charac- 
terized by q g and q s ), but there is a discontinuity at 
the transition size D t with the strength scaled compo- 
nent shifted down by an appropriate amount x t (see their 
Fig. 3b). This means that / max should be higher than 
that derived using equation (16) with q = q g by a factor 



1) of initial mass M to t(0) = lM e 
= 2700 kg m" 3 , e = 0.05, e/I = 1, 



x t (3( ? g-5)(3(Z s -5)- 1 pbi/A) 3(9 ^ 9s) . Since x t < 1, then 
substituting the values from Benz & Asphaug (1999) 
given above implies that Table 1 underestimates / max by 
at most a factor of 50-100 (possibly much less). In other 
words, we anticipate that by including a more realistic 
prescription for and the resulting size distribution, 
this would change the inferred / ma x but not upwards by 
an amount more than two orders of magnitude. For this 
reason, transience is only inferred for those systems for 

Which / b s / j 'max > 100. 

A lower eccentricity is, however, one potential av- 
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enue for increasing the amount of dust remaining at late 
times. Equation (20) shows that since G(ll/6, X c ) oc 
e 5/3 (pjg ^a), this means that reducing e from 0.05 
to 0.01 or 0.001 gives a decrease in G(U/6,X C ) of 15 
or 680, and so an increase in / max by a corresponding 
amount (e.g., Fig. 3). In fact the increase can be much 
more than this, since when e is reduced to levels below 



tant when v IC \ < v c 



v /(2/3)7r / 9G(10- 3 D) and so 



4.7xlO- 5 y / g*rM*- 1 [1.25+ (I/e) 2 }- 1 then A c > 1 (e.g., 

Fig. lb). In such a regime mutual collisions do not result 
in the destruction of planetesimals, rather in their merger 
and growth. At this point G(ll/6, X c ) < 0, i.e., / max is 
infinite since, in this simple model, whatever the start- 
ing conditions there is no evolution (although in practice 
the size distribution would evolve due to planetesimal 
growth). At <~ 1 AU, this means e must be larger than 
0.0005 (for Q* = 200 J kg -1 , appropriate for D c = 0.15 
km) or 0.014 (for = 2 x 10 5 J kg -1 , appropriate for 
D c = 2000 km) to initiate a collisional cascade, values 
which are consistent with those quoted by more detailed 
planet formation models (e.g., Kenyon & Bromley 2002). 
Such eccentricities would be expected through stirring 
either by > 1000 km planetesimals which formed within 
the belt, or by more massive perturbers which formed 
outside the belt, both of which can be expected to occur 
within 10-100 Myr (Kenyon & Bromley 2006). This was 
considered in §3.1 where it was shown that the cascade 
would be initiated following the growth of ~ 2000 km 
planetesimals on timescales that are much shorter than 
the age of the system for the disk masses required to 
produce a dust luminosity at the observed level. 

The only route which could plausibly maintain the hot 
dust systems in Table 1 in collisional equilibrium over 
the age of the stars might be to invoke some mecha- 
nism which maintains the eccentricity at a level which 
the cascade is only just being eroded. However, Fig. la 
shows that G(11/6,A C ) is a strong function of e when 
G(11/6,X C ) < 1, since the range G(11/6,A C ) = 0-1 
is covered by a factor of less than two in eccentricity. 
Thus we consider it reasonable to assume that the best 
possible combination of and e in this regard would 
result in G(ll/6, X c ) w 1 (corresponding to X c = 0.69); 
lower values of G(ll/6, X c ) are possible, but only within 
a very narrow range of eccentricity. Since in the above 
example with a realistic Qq prescription extending up to 
2000 km we assumed e = 0.05 which already resulted in 
G(q g ,X c ) < 1, we consider that it is not reasonable to 
fine tune the eccentricity further to increase / max ; e.g., 
decreasing to e = 0.03 results in some disks not evolving 
and the rest with / max higher than that quoted in Table 
1 by a factor ~ 150. Thus we conclude that the estimate 
given in Table 1 (and e.g., equation 16) underestimates 
/max by at most a factor of ~ 100, unless the eccentricity 
happens to lie within ±10% of a critical value. 

It is also worth noting that low levels of eccentricity 
would result in large gravitational focussing factors for 
large planetesimals which would enhance / cc and so de- 
crease the time for these planetesimals to be catastroph- 
ically destroyed, something which is compounded by the 
higher collision velocity in gravitationally focussed colli- 
sions which reduces X c because collisions with smaller 
planetesimals can cause catastrophic disruption (e.g., 
equation 11). However, we do not need to account for 
this here, since gravitational focussing becomes impor- 



when e < 4 x 10~ 7 ^ prM* 1 [1.25 + (I/e) 2 ]~ 1 D (where 

D is in km); i.e., when e < 2 x 10~ 6 for D = 0.15 km 
and e < 0.027 for D = 2000 km at 1 AU from a 1M Q star, 
both of which occur close to or below the level at which 
collisions result in accumulation rather than destruction. 

3.3. Are these the products of single collisions? 

One possible origin for the hot dust which is quoted in 
the literature is that it is the product of a single colli- 
sion (Song et al. 2005). Our model can be used to make 
further predictions for the likelihood of massive collisions 
occurring within an asteroid belt. The maximum num- 
ber of parent bodies (i.e., planetesimals) larger than D pb 



remaining at late times occurs when M to t 
so is given by: 



Mr, 



and 



n(D>D ph )- 



5.6 x 10 w r 3 - 5 (dr/r) 

m - 5 n 2 t 

(3<7-5)[l + 1.25(e/7) 



21-0.5 



3<j-3 



(3-Sq)G(q,X c ) 

The collision timescale for planetesimals of size D p b is 

t c (D pb )=t c (D c )f cc (D c )/f cc (D ph ) 

= 10 6 t agc (D ph /D c ) 3 «- 5 , (23) 

noting that the collisional lifetime of the largest plan- 
etesimals, t c (D c ) 1 is the age of the star for a planetes- 
imal belt at maximum luminosity for this age. These 
can be combined to give the destructive collision rate for 
planetesimals larger than D^: 



(22) 



dN c {D > D ph )/dt = l000r 13 / 3 {dr/r)t- 2 D c D. 



P b 



M^ 4/3 Q^V 5 / 3 , (24) 

in Myr -1 , where the assumptions that g = 11/6, e = I 
and X c < 1 have been used in deriving this equation. 

We now assume that we are considering collisions capa- 
ble of reproducing the observed dust level, / bs, so that 
the lifetime of the resulting collision products can be esti- 
mated from the collisional lifetime of that dust, assumed 
to be of size D u (WD02): 



t c (D u ) = 0.04r i5 M- u ' 5 (dr/r)/ o - 



bs> 



(25) 

in years, noting that collisions would remove the dust on 
a faster timescale than P-R drag (Wyatt 2005; Beichman 
et al. 2005). Combining equations (24) and (25) gives 
the fraction of time that collisions are expected to result 
in dust above a given level of /obs^ 



P(f > Jobs) =4 x \Q-\ 3 ^(dr/rft- 2 c D c D- 3 x 



11/6 



fobsQii 



5/6-5/3 



(26) 



To estimate the minimum size of the parent body, D pri , 
responsible for this dust, we consider how large a plan- 
etesimal must be to reproduce / Q b s if a destructive colli- 
sion resulted in one fragment with half the mass of the 
original planetesimal (i.e., the definition of a destructive 
collision), with the remaining mass in particles of size 
Du- 

D ph = 890[£ bl r 2 /o bs ] 1/3 . (27) 



10 



Wyatt et al. 



TABLE 2 

Parameters in the model for the hot dust systems of Table 1 used to determine whether the observed dust can be the 
outcome of a single collision in a massive asteroid belt which is itself not normally bright enough to be detected. 



Star name 


D pb , km 


N(D > D pb ) 


dN c (D > D ph )/dt, Myr- 1 


tc(D hl ), yr 


P(f > /obs) 


Single collision? 


HD98800 


530 


200 


41 


0.36 


15 x lO" 6 


No 


HD113766 


280 


890 


150 


41 


6100 x 10~ 6 


Not imposs 


HD12039 


110 


77,000 


12,000 


2300 


27* 


Not imposs 


BD+20307 


320 


0.47* 


0.0039 


0.49 


0.0019 x 10~ 6 


No 


HD72905 


15 


1.4 


0.036 


22 


0.79 x 10~ 6 


No 


r\ Corvi 


110 


7.8 


0.033 


59 


2.0 x 10~ 6 


No 


HD69830 


39 


19 


0.068 


110 


7.7 x 10~ 6 


No 



For disks with P(f > / b s ) > 1, this value indicates the number of collisions at that level we can expect to sec in the disk at any one time. 
Likewise, for disks with N(D > -Dpb) < 1, this value indicates the probability that there is an object of this size remaining in the disk. 



Table 2 lists the parameters for the hot dust systems 
assuming the canonical parameters of = 200 J kg -1 , 
D c = 2000 km and e = 0.05. To determine whether 
a system could have been reproduced by a single colli- 
sion, the final value of P(f > /obs) was compared with 
the statistic that 2% of systems exhibit hot dust (which 
therefore considers the optimistic case where all stars 
have planetesimal belts at a few AU). For the systems 
which were inferred in Table 1 to be transient, all arc 
extremely unlikely (< 0.001%) to have been caused by a 
single collision amongst planetesimals in a planetesimal 
belt which has undergone a collisional cascade since the 
star was born. 

While this statistic is subject to the uncertainties in 
the model parameters described in §3.2, and so could 
be in error by around two orders of magnitude, it must 
also be remembered that the most optimistic assump- 
tions were used to arrive at this figure. For example, it 
is unlikely that the destruction of planetesimals of size 
Z?pb would release half of the mass of the planetesimal 
into dust D^i in size. 5 On the other hand, one might 
consider that the lifetime of the observed dust, i c (Z?bi), 
is an underestimate of the duration of dust at the level 
°f / > /obs, since the dust could be replenished from 
the destruction of larger particles. Indeed Farley et al. 
(2006) model the destruction of a 150 km planetesimal 
in the asteroid belt and inferred a dust peak that lasted 
~ 1 Myr, precisely because large fragments produced in 
the collision replenished the dust population. However, 
it should be cautioned that the dust peak inferred by 
Farley et al. (2006) would not have been detectable as 
an infrared excess since it only caused a factor ~ 10 en- 
hancement in the luminosity of the zodiacal cloud (i.e., 
to / w 0.8 x 10~ 6 ), and that in the context of our model, 
invoking a population of larger grains that result from 
the collision would lead to a larger parent body (i.e., a 
larger D p b) required to reproduce the observed luminos- 
ity /obs and so less frequent collisions; i.e., it may be 
possible (even desirable) to increase t c (D D i), but only 
at the expense of decreasing dN c (D > D p b)/dt leading 
to little change in P(f > /obs)- We note that £ C (-Dbi) 

5 Such an optimistic assumption should not be dismissed out of 
hand, however, since the large amount of collisional processing that 
must have taken place means that planetesimals more than a few 
km would be rubble piles. These would have undergone shattering 
and reaccumulation numerous times meaning that they could have 
deep dusty regolith layers which could be preferentially ejected in 
a collision. 



given in Table (2) is sufficiently short that a measure- 
ment of the variability of the infrared excess on realistic 
(few year) timescales could lead to constraints on the size 
of the grains feeding the observed phenomenon, since if 
a population of larger grains existed then the luminosity 
would fade on much longer timescales. 

A further argument against the transient disks being 
caused by single collisions is the fact that the probability 
of seeing the outcome of a collision, P(f > /obs), falls off 
oc t~g C , which means that we would expect to see more 
transient disks around younger stars than around older 
stars (because young stars have more massive disks with 
more large planetesimals and so more frequent collisions) . 
There is some evidence from Table 1 that transience is 
more common around young systems, since none of the 
transient systems is older than 2 Gyr, whereas sun-like 
stars in the solar neighborhood would be expected to 
have a mean age of ~ 5 Gyr. However, while the statistics 
are poor, a t~^ e dependence does seem to be ruled out; 
e.g., we would have expected to have detected 10 times 
more transient disks caused by single collisions in the 
age range 50-500 Myr 6 than in the age range 0.5-5 Gyr, 
whereas 2 transient disks are known in the younger age 
bin, and 2 in the older age bin which is more consistent 
with a t~g e dependence. 

In fact, within the context of this model, all of the disks 
which we infer to be transient would also be inferred to 
not be the product of single collisions. This is evident 
by substituting D p b from equation (27) and / max from 
equation (18) into equation (26) to get: 

P(f > /obs) = 0.2 x lOVnW/obs^M^r-^- 1 ) 5 / 6 , 

(28) 

which reduces to P(f > /obs) = 
16(/max//obs) 2 (A^*f ~ 1 ) 5//6 for the canonical param- 
eters used before. Since transient disks are defined 
by /obs//max ^> 1000, this means they cannot also 
have a high probability of having their origin in single 
collisions. It would only be inferred that disks with 
/obs//max *C 1000 could have their origin in single 
collisions, but since it is also possible that these disks 
are the result of steady state collisional evolution there 
is no need to invoke a single collision to explain their 
presence, which is why Table 2 simply concluded that 
it is "not impossible" that the disks of HD113766 and 

6 It is not reasonable to extend the age range to younger systems, 
since, as noted in §3.1 it is hard to discern whether or not dust 
detected in such systems is transient. 
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TABLE 3 

Parameters in the model for the transient hot dust systems of Table 1 used to determine whether the observed dust 
could originate in the destruction of a planetesimal belt coincident with the dust: dm loss /dt is the observed mass loss 
rate, m max is the maximum mass of a planetesimal belt that is coincident with the dust given the age of the star, 

*(/ > fobs) IS THE LENGTH OF TIME SUCH A PLANETESIMAL BELT COULD SUSTAIN THE OBSERVED DUST LUMINOSITY, AND r out (100Myr) IS 

the radius of a planetesimal belt which would still have enough mass to sustain the observed dust luminosity for 100 

Myr. 



Star name 


dM loss /dt, lO^MeMyr" 1 


M m ax, 10- 6 M ffi 


t(.f > fobs), Myr 


rout(lOOMyr), AU 


BD+20307 


8.0 X 10 6 


53 


6.7 X 10~ 6 


45 


HD72905 


19 


0.072 


3.7 x 10~ 3 


2.4 


r) Corvi 


2500 


57 


0.023 


9.6 


HD69830 


64 


12 


0.18 


4.5 



HD12039 are the product of single collisions. What 
equation (28) does indicate, however, is that it is 
possible for single collisions to cause disks to spend some 
fraction of their time at a luminosity enhanced above 
the nominal maximum value / max , and that this occurs 
more readily for disks at smaller radii and around higher 
mass stars. However, whether single collisions really do 
achieve an observable increase in luminosity depends 
on the size distribution of the collisional fragments, for 
which it must be remembered that equation (28) used 
an unrealistically optimistic estimate. 

3.4. Are parent planetesimals coincident with dust? 

For similar reasons to those in §3.3 it is also possible 
to show that the parent planetesimals of the dust are 
extremely unlikely to originate in a planetesimal belt that 
is coincident with the dust. The reason is that the mass 
remaining in such a belt would be insufficient to replenish 
the dust for a length of time commensurate with the 
statistic that 2% of stars show this phenomenon. The 
observed dust luminosity, assuming this is comprised of 
dust of size -D D i which has a lifetime of t c (Db,\) (equation 
25), implies a mass loss rate due to mutual collisions 
between the dust grains of: 

dM loss /dt = 1700/ o 2 bs r°- 5 ^M-°- 5 (r/dr), (29) 

in M^/Myr, and this is independent of the collisional 
evolution model of §2. However, due to the collisional 
evolution of a planetesimal belt's largest members, there 
is a maximum mass that can remain in a belt at the 
same radius as the dust at this age, and this is given in 
equation (15). This means that if the observed dust orig- 
inates in an event which, for whatever reason, is causing 
planetesimals in a belt at the same radius as the dust 
to be converted into dust, then this can last a maximum 
time of t(f > /obs) = M max / dM\ oss /dt before the plan- 
etesimal belt is completely exhausted. These figures are 
given in Table 3 which shows that the longest the type 
of transient event observed could be sustained in these 
systems is under 1 Myr, under the assumptions about 
the planetesimal belts employed in the rest of the paper. 

A maximum duration of 1 Myr is not sufficient to ex- 
plain the statistic that 2% of sun-like stars exhibit this 
phenomenon, since the median age of such stars is 5 
Gyr, indicating a typical duration (even if this occurs 
in multiple, shorter, events) of around 100 Myr. Clearly 
a reservoir of mass is required in excess of that which it 
is possible to retain so close to the star. 

3.5. Constraints on parent planetesimal belt 



If we assume that the observed mass of hot dust orig- 
inates in planetesimals that were initially in a belt at a 
radius r out which has properties like those assumed in 
the rest of the paper, and a fractional luminosity of / ou t, 
then there are two main constraints on that belt. First, 
assuming that this belt has been collisionally evolving for 
the age of the star, then this belt cannot have more mass 
(or luminosity) than the maximum that could possibly 
remain due to collisional processing, i.e., / ou t < /max 
(equation 16). Second, it must have sufficient mass re- 
maining to feed the observed mass loss rate for long 
enough to reproduce the statistic that 2% of stars ex- 
hibit this phenomenon which implies a total duration of 
> 100 Myr. For a belt to have enough mass to feed the 
observed hot dust luminosity of / D b s at a radius r for 
a total time of thot m Myr requires the belt to have a 
luminosity of: 

/out > 710t hot /2 bB r-V°- B I> c - - 5 i2- B (dr/r)- 1 , (30) 

or rather, this is the luminosity it must have had before 
it was depleted. 

Comparing this with the maximum mass possible at 
this age indicates that the parent belt must have a min- 
imum radius of: 

rout(thot) >615^ t 13 ^ c 13 / 6 b / s 13 r 3/26 (dr/r)- 6 / 13 x 

D -3/13g*-5/26 e 5/13 £ 3/13 M B/26_ ( g 1} 

Table 3 gives an estimate of the minimum radial location 
of such a planetesimal belt, under the assumption that 
the event (or multiple events) of high hot dust luminosity 
last ihot = 100 Myr. These values indicate that the plan- 
etesimal belts must be at least a few AU from the star. 
It must be cautioned that this conclusion is relatively 
weak in the case of HD69830, since the uncertainty in 
the properties of the planetesimals still leaves two orders 
of magnitude uncertainty in the maximum luminosity, 
/max, and so also in the maximum mass M max (see §3.2). 
This means that, with suitable planetesimal belt prop- 
erties, a belt in this system that is coincident with the 
dust at 1 AU may be able to replenish the observed phe- 
nomenon for 20 Myr. However, we still consider this to 
be an unlikely scenario, since it would require that the 
mass of the planetesimal belt is depleted at a constant 
rate for the full 100 Myr, whereas most conceivable sce- 
narios would result in a mass loss rate which decreases 
with time as the planetesimal population is depleted thus 
requiring an even larger starting mass. 

These two constraints are summarized for the 4 sys- 
tems with transient hot dust in Fig. 4, which shows the 



12 



Wyatt et al. 





0"' 



cr'g 



i HD72905 


f m0) / 


: X. ^\ 70u.m detection 
: V ^^irrfit (BB) 




- O N. ^\ 

i hot dust N. \^ 


cold dust / 


/ moss loss limit^^ 


possible porent 
planetesimal 
belt 




Fig. 4. — Constraints on the fractional luminosity and radius of the planetesimal belt feeding the observed transient hot dust (shaded 
region) for the systems: (top left) r\ Corvi, (top right) HD72905, (bottom left) HD69830 and (bottom right) BD+20307. The solid 
lines are the constraints imposed by the far-IR detection limits (assuming black body emission), the maximum luminosity possible in a belt 
at this radius due to erosion by collisional processing, and the luminosity from a belt of sufficient mass to feed the observed mass loss rate 
for 100 Myr. The properties of the hot dust in these systems is shown with a diamond, and those of the cold dust, where known (the top 
two figures), shown with a triangle. 



shaded region of parameter space in / out and r out where 
the parent planetesimal belt can lie. This figure also 
shows the location of the hot dust at / Q b s and r, illus- 
trating the conclusion of §3.1 that this lies significantly 
above f max , the maximum fractional luminosity expected 
for a planetesimal belt at the age of the parent star. Note 
that the value r out (100Myr) given in Table 3 denotes the 
intersection of the limits from / max and from equation 
(30). 

A third constraint for the parent planetesimal belt 
comes from far-IR observations of these systems. For 
2/4 of the transient dust systems a colder dust compo- 
nent has already been detected: rj Corvi has a planetes- 
imal belt with a resolved radius of ~ 100 AU (Wyatt et 
al. 2005), and HD72905 has one inferred to be at - 14 
AU (Beichman et al. 2006a). In both cases these outer 
planetesimal belts have been inferred to be at a differ- 
ent spatial location from the hot dust either because of 
imaging constraints (Wyatt et al. 2005) or from anal- 
ysis of the SED (Beichman et al. 2006a). The proper- 
ties inferred for these planetesimal belts are indicated on 
Fig. 4 and lie within the shaded region, implying that 
these planetesimal belts do not have to be transiently 
regenerated, and also provide a plausible source popula- 
tion for the hot dust found closer in. However, no such 
excess emission has been seen toward HD69830 at either 
70 ^m (Beichman et al. 2005) or 850 fim (Sherct, Dent & 



Wyatt 2004) indicating a planetesimal belt with a mass 
at most 5-50 times greater than our own Kuiper belt. 
Likewise, BD+20307 does not have a detectable excess 
in IRAS 60 ^m observations (Song et al. 2005). 

A low mass reservoir of planetesimals does not neces- 
sarily rule out the presence of an outer planetesimal belt 
which is feeding the hot dust for two reasons. First, the 
shaded region of Fig. 4 actually constrains the proper- 
ties of the planetesimal belt at the time at which de- 
pletion started; i.e., this population may have already 
been severely depleted by the same event which is pro- 
ducing the dust and we are now nearing the end of the 
hot dust episode. Second, the constraints imposed by a 
non-detection in the far-IR do not eliminate the whole 
of the parameter space in which an outer planetesimal 
belt can lie. Fig. 4 includes the constraints on the outer 
planetesimal belt imposed by the non-detection of excess 
in the far-IR, assuming that the dust emits like a black 
body. The resulting detection limit is then given by: 

/det = 3.4 x 10 9 F dct (A)d 2 r ou 2 t /B 1/ (A, T bb ), (32) 

where F^et is the detection limit in Jy d is the distance to 
the star in pc, and B v {\, T bb ) is in Jy/sr. For BD+20307 
the non-detection is limited by the sensitivity of IRAS, 
and so lower limits should be achievable with Spitzer. For 
the two systems with non-detections, the shaded region 
already takes the far-IR constraint into account. 
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The simplification that the emission comes from black 
body-type grains means that equation (32) underesti- 
mates the upper limit from the far-IR fluxes. This is 
because the majority of the luminosity comes from small 
grains which emit inefficiently at long wavelengths. In- 
deed, the black body assumption would require the hot 
dust of HD69830 and BD+20307 to have been detected 
in the far-IR, whereas this is not the case. We modeled 
the emission from non-porous silicate-organic refractory 
grains in a collisional cascade size distribution at 1 AU 
from these stars to find that the black body assumption 
used in equation (32) underestimates the limit by a fac- 
tor of 3 — 5 meaning that non-detection of the hot dust 
in these systems in the far-IR is to be expected. This 
also means that slightly more luminous outer planetesi- 
mal belts than those indicated by the shaded region on 
Fig. 4 may still have escaped detection in the far-IR. 

Until now we have not proposed a mechanism which 
converts the planetesimals into dust. Whereas Beichman 
et al. (2005) invoke sublimation of comets as the origin 
of the hot dust, and use this to estimate the mass of the 
parent planetesimal belt, we consider a scenario in which 
a significant fraction of material of all sizes in the parent 
planetesimal belt is placed on orbits either entirely coin- 
cident with the hot dust, or with pericenters at that loca- 
tion. In this scenario the dust is reproduced in collisions 
and the material maintains a collisional cascade size dis- 
tribution. Simply moving material from r out to r would 
result in an increase in fractional luminosity from / out to 
/out(foutA) 2 • This indicates that the parent planetesi- 
mal belt responsible for the hot dust could have originally 
been on the line on Fig. 4 traced by / ou t = /obs^Aout) 2 - 
Since this is parallel to the mass loss limit line (equa- 
tion 30), and for all but HD69830 the observed hot dust 
component lies below this line, this indicates that par- 
ent planetesimal belts in the shaded region could be re- 
sponsible for the hot dust observed, as long as a large 
fraction of their mass is scattered in to the inner regions. 
However, it is to be expected that only a fraction of the 
outer planetesimal belt ends up in the hot dust region, 
and so it is more likely that the parent planetesimal belt 
started on a line which falls off less steeply than oc r m 2 t , 
and this is consistent with the ratio of the hot and cold 
components of r\ Corvi and HD72905 which indicate a 
dependence of / out = /obs(r/r out )°- 5±0 - 2 ; it is also inter- 
esting to note that both have r out /r = 60 — 70. We defer 
further consideration of the expected properties of the 
parent planetesimal belt to a more detailed model of the 
dynamics of the types of events which could cause such a 
perturbation, but simply note here that the existence of 
an outer planetesimal belt is not ruled out by the current 
observational constraints in any of the systems. 

4. DISCUSSION 

A simple model for the steady state evolution of dust 
luminosity for planetesimal belts evolving due to colli- 
sions was described in §2. This showed how at late times 
the remaining planetesimal belt mass and so dust lumi- 
nosity is independent of the initial mass of the belt. This 
has important implications for the interpretation of the 
properties of detected disks. This paper discussed the 
implications for the population of sun-like stars with hot 
dust at < 10 AU; the implications for the statistics will 
be discussed in a forthcoming paper (Wyatt et al., in 



prep.). 

It was shown in §3.1 that for 4/7 of the systems with 
hot dust their radius and age are incompatible with a 
planetesimal belt which has been evolving in quasi-steady 
state over the full age of the star, and in §3.2 it was 
shown that this is the case even when uncertainties in 
the model are taken into account. This implies either 
that the cascade was started recently (within the last 
Myr or so), or that the dust arises from some other tran- 
sient event. Recent ignition of the collisional cascade 
seems unlikely, since the mass required to feed the ob- 
served luminosity would result in the growth of 2000 km 
planetesimals which would stir the belt and ignite the 
cascade on timescales much shorter than the age of the 
stars. Possible origins for the transient event that have 
been proposed in the literature are: recent collision be- 
tween massive planetesimals in a planetesimal belt which 
introduces dust with a size distribution <; > 11/6 and so 
can be detected above a collisional cascade which is too 
faint to detect; one supercomet ~ 2000 km in diameter 
that was captured into a circular orbit in the inner sys- 
tem replenishing the dust through sublimation (Beich- 
man et al. 2005); a swarm of comets scattered in from 
the outer reaches of the system (Beichman et al. 2005). 
In §3.3 the collisional model was used to show that the 
transient disks are very unlikely (< 0.001% for the most 
optimistic estimate for any of the stars compared with 
a detection probability of 2% for transient hot dust) to 
have their origin in a recent collision; such collisions oc- 
cur too infrequently. In §3.4 it was also shown that the 
parent planetesimals of the observed dust must originate 
in a planetesimal belt much further from the star than 
the observed dust, typically at 3> 2 AU. This is because 
collisional processing means that the mass that can re- 
main so close to the star at late times is insufficient to 
feed the observed phenomenon. 

The most likely scenario is thus a recent event which 
provoked one or more planetesimals to be scattered in 
from further out in the disk (Beichman et al. 2005). 
The observed dust could have been produced from such 
a scattered planetesimal population through their grind- 
ing down in mutual collisions (§3.5), although sublima- 
tion close to the pericenters of the planetesimals' orbits 
is a further possible source of dust. More detailed study 
of the scattering and consequent dust production pro- 
cesses is required to assess these possibilities. However, 
this scenario is supported by the presence of far-IR emis- 
sion originating from a colder outer planetesimal belt 
component in 2/4 of the transient dust systems. The 
constraints on the outer planetesimal belt which is feed- 
ing the phenomenon are discussed in §3.5, showing that 
the outer planetesimal belts already found in i] Corvi 
and HD72905 provide a plausible source population for 
the hot dust found closer in, and that the current non- 
detection of cold dust around the remaining two systems 
does not rule out the presence of an outer planetesimal 
belt capable of feeding the observed hot dust luminosity. 

One clue to the origin of the parent planetesimals of 
the dust may be the composition of that dust. Silicate 
features have been detected in the mid-IR spectrum of 
all of the transient hot dust stars (Song et al. 2005; Be- 
ichman et al. 2005; Beichman et al. 2006a; Chen et al. 
2006). Detailed modeling of the spectrum of HD69830 
indicates that the mineralogical composition of its dust is 



14 



Wyatt et al. 



substantially different from that of comets, rather there 
is a close match to the composition of P- or D-type as- 
teroids found mainly in the 3-5 AU region of the solar 
system (Lisse et al. 2006). While the radial location 
at which planctesimals of this composition form in the 
HD69830 system will depend on the properties of its pro- 
tostellar nebula, which may be significantly different to 
that of the protosolar nebula, as well as on the structure 
and evolution of its planetary system, evidence for water 
ice in the dust spectrum indicates that the parent body 
formed beyond the ice-line in this system (Lisse et al. 
2006), i.e., beyond 2- 5.5 AU (Lecar et al. 2006; Alibert 
et al. 2006). Thus the compositional data supports the 
conclusion that the dust is not produced by a planetesi- 
mal that formed in situ. However, it is worth noting that 
the same compositional data also finds evidence for dif- 
ferentiation in the parent body (inferred from abundance 
differences between the dust and the star) and for heat- 
ing of its rocky material to > 900 K (inferred from the 
absence of amorphous pyroxene), which would also have 
to be explained in the context of an outer planetesimal 
belt origin for the dust. 

An analogous transient event is thought to have hap- 
pened in the solar system resulting in the period known 
as the Late Heavy Bombardment (LHB) when the terres- 
trial planets were subjected to an abnormally high im- 
pact rate from asteroids and comets. This is believed 
to have been triggered by a dynamical instability in 
the planetary system resulting from Jupiter and Saturn 
crossing the 1:2 resonance during their slow migration 
(inwards for Jupiter, outwards for Saturn) due to angu- 
lar momentum exchange with the primordial Kuiper belt 
(Gomes et al. 2005). In this scenario both the asteroid 
and Kuiper belts were depleted with a large fraction of 
these objects being scattered into the terrestrial planet 
region during an event which lasted 10—150 Myr (Gomes 
et al. 2005), i.e., exactly the type of event required to 
explain the observed hot dust in the scenario proposed 
here (§3.5). Dynamical instabilities in cxtrasolar plan- 
etary systems can also arise from mutual gravitational 
perturbations between giant planets which formed close 
together (Lin & Ida 1997; Thommes, Duncan & Levison 
1999). In both scenarios slow diffusion of the orbits of 
the planets means that the dynamical instability can oc- 
cur up to several Gyr after the formation of the planetary 
system. The delay to the onset of the instability is de- 
termined by the separation of the outer planet from the 
outer planetesimal belt (Gomes et al. 2005), or from the 
separation between the planets (Lin & Ida 1997), with 
larger separations resulting in longer timescales. 

Little is known about the planetary systems of four of 
the hot dust systems. However, three Neptune mass (or 
Jupiter mass if the system is seen face-on) planets have 
recently been discovered orbiting the star HD69830 at 
< 1 AU on nearly circular orbits (Lovis et al. 2006). 
Dynamical simulations showed that the detected plan- 
etary system is stable on timescales of 1 Gyr. This 
does not, however, rule out the possibility of a dynam- 
ical instability having occurred. While no mean mo- 
tion or secular resonances are immediately identifiable 
within the detected planetary system which could have 
have been crossed recently invoking such a catastrophic 
event, it is possible that the instability arose with an- 



other planet further out which has yet to be detected 
with longer timescale observations. It is also possible 
that a fourth planet which existed in the region 0.19- 
0.63 AU between the planets HD69830c and HD69830d 
has recently been scattered out due to a dynamical insta- 
bility (e.g., Thommes et al. 1999). The region 0.3-0.5 AU 
was identified in Lovis et al. (2006) as being marginally 
stable, and to encompass several mean motion resonances 
with the outer planet, including the 1:2 resonance at 0.4 
AU; i.e., a putative fourth planet could have remained 
in this region for the past 2 Gyr until the slow migra- 
tion/diffusion of the outer planet (HD69830d) caused the 
1:2 resonance to coincide with the orbit of the putative 
planet which was then scattered outward thus promot- 
ing the depletion of an outer planetesimal belt much of 
which was scattered into the inner regions of the system. 
Alibert et al. (2006) considered that the most plausible 
formation scenario for the planetary system of HD69830 
included the inward migration of the outer planets from 
beyond the ice-line at a few AU. This would put a sub- 
stantial distance between the outer planet (HD69830d) 
and any outer planetesimal belt which favors a delay of 
2 Gyr before the onset of the instability. Searches for 
further planetary companions in this system, and for the 
relic of its outer planetesimal belt, are clearly necessary 
to constrain the evolutionary history of this system. 

In conclusion, ~ 2% of sun-like stars exhibit transient 
hot dust in the terrestrial planet region; this dust must 
originate in a planetesimal belt located further from the 
star than the dust, typically at » 2 AU. Just four mem- 
bers of this class are currently known, although it seems 
reasonable to assume that our own solar system would 
have been placed in this class during the LHB. The fre- 
quency of this phenomenon indicates that either all stars 
are subjected to an epoch of similar duration (lasting 
~ 100 Myr assuming a typical age of 5 Gyr) or that a 
smaller fraction of stars undergo much longer (or multi- 
ple) events. The distribution of the ages of the stars in 
this class indicate that the likelihood of these events oc- 
curring falls off roughly inversely proportional to the age 
of the stars. An origin for these events in a dynamical 
instability as proposed for the LHB in the solar system 
is supported by the recent discovery of a multiple planet 
system coincident with the dust in one of the systems 
currently in this class. However, since the LHB in the 
solar system is thought to have lasted just ~ 100 Myr, it 
remains to be seen whether we are to infer that dynami- 
cally unstable planetary systems form around all stars, or 
that the LHB event in other systems lasted much longer 
than in our own, or perhaps that there is in fact more 
than one mechanism causing this hot dust signature. Ob- 
servations that further constrain the planet, planetesimal 
and dust complements of the transient hot dust systems 
are needed to ascertain the similarities and dissimilarities 
within this population. 
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TABLE Al 

Symbols employed in this paper and their units. 



Symbol 


Units 


Meaning 


B V (\,T) 


Jy/sr 


black body emission spectrum 


d 


pc 


distance to star 


dM loss /dt 


M e /Myr 


rate of mass loss assuming observed dust has size D D i 


dr 


AU 


planetesimal belt width 


D M 


(jm 


diameter of smallest dust in cascade 


D c 


km 


diameter of largest planetesimal in cascade 


D cc 


km 


smallest planetesimal capable of destroying planctcsimals of diameter D c 


D pb 


km 


minimum diameter of parent body required to produce observed dust 


D t 


km 


planetesimal diameter at transition between strength and gravity regimes 


e 




mean orbital eccentricity of planetesimals 


f 




fractional luminosity (= Li r /L*) in model 


/det 




fractional luminosity for emission from belt to be detected 


/max 




maximum fractional luminosity of cascade after time t age 


Jobs 




fractional luminosity observed 


/out 




fractional luminosity of putative outer planetesimal belt feeding the dust 


f(e,I) 




ratio of collision velocity to Keplerian velocity 


fee 




see equation (8) 


FdetW 


Jy 


detection limit at wavelength A 


G(q,X c ) 




sec equation (9) 


I 


rad 


mean orbital inclination of planctcsimals 


K 




scaling factor in size distribution 




L e 


stellar luminosity 




Lq 


' c 1 1 ■ ■ i e j • i • 1 1 i 

mtrarcd luminosity ot material m the cascade 




Mq 

M® 


stellar mass 


M ma i 


maximum mass remaining in cascade after time t age 


Mtot 


M© 


total mass of material in cascade 


n(D) 




size distribution ot material in the cascade 


n(D > Dpb) 




1 C 1 • J_ • 11 J_l T~\ 

number ol objects in cascade larger than D va 


dN c (D > D ph )/dt 


IT 1 

Myr 1 


destructive collision rate for planctcsimals larger than Dpi-, 


P(f > fobs) 




fraction of time collisions result in / > f \, s 


1 

Is 




slope of size distribution 

slope of size distribution expected in the gravity regime 


q s 




slope of size distribution expected m the strength regime 


Qi> 


J kg 


specific incident energy required to catastrophically destroy a pfanctcsimal 


r 


ATT 

AU 


planctcsimaf beft radius, assumed to be coincident with dust 


rout {that ) 


AU 


outer planctcsimaf beft radius required to maintain / b s for r no t 


Tout 


ATT 

AU 


radius of putative outer planetesimal belt feeding the dust 






cxponciiL in iciaitioii "t/j-j u 


t(f > fr.U B ) 


Myr 


time a planetesimal belt at v can sustain f ^> fobs 


iage 


Myr 


time since cascade initiated (assumed to be stellar age) 


*c 
tc(D) 


yr 


collisional lifetime of planctcsimals of size D c 


yr 


collisional lifetime of material of size D 


that 


Myr 


total duration of hot episodes throughout stellar lifetime 


Vre\ 


m s _1 


relative velocity of collisions 


Vcsc 


m s _1 


escape velocity 


Vk 


m s _1 


Keplerian velocity 


Xt 




jump in size distribution expected at Dt 


x c 




= D cc /D c 


P 


kg m -3 


planetesimal density 




AU 2 


total cross-sectional area of material in cascade 



APPENDIX 

SUMMARY OF SYMBOLS 

The symbols which are employed in this paper are summarized in Table Al along with the units assumed throughout 
the paper. 
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